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I Abstract 

' Search of a proper and reahstic equations of state (EOS) for 

^ ' strongly interacting matter used in the study of QCD phase diagram 

, still appears as a challenging task. Recently, we have constructed a 

' hybrid model description for the quark gluon plasma (QGP) as well as 

hadron gas (HG) phases where we use a new excluded-volume model 
for HG and a thermodynamically-consistent quasiparticle model for 
the QGP phase. We attempt to use them to get a QCD phase bound- 
, ary and a critical point. We test our hybrid model by reproducing 

the entire lattice QCD data for strongly interacting matter at zero 
baryon chemical potential (/j,s)and predict the results at finite hb 
and T. 



1 Introduction 



Quantum chromodynamics (QCD) predicts that at sufficiently high tem- 
Q>^ ■ peratures (T) and/or chemical potentials (/is), strongly interacting matter 

I goes through a phase transition from colour insulating hadron gas (HG) 

phase to colour conducting quark gluon plasma (QGP) phase [l][2]. Heavy 
! ion collisions can provide a unique opportunity to study this QCD phase 

■ transition from HG to QGP. Our endeavour in such studies is the search 
I of a suitable equation of state (EOS) for the description of both phases of 

strongly interacting matter. Significant success has been gained in lattice 
calculations using QCD thermodynamics to provide a valid EOS for QCD 
matter at zero baryon chemical potential but at large T. However, the 
lattice methods are unreliable to describe the properties of matter at finite 

■ density of baryons. Therefore, finding an EOS for QCD matter valid at 
zero as well as non-zero chemical potential is a challenging problem. In this 
paper, we construct a hybrid model where a new excluded-volume model 
is proposed for HG description and a thermodynamically-consistent quasi- 
particle model is given for the QGP phase in predicting the properties of 
entire QCD matter. Moreover, we use the hybrid model to get the QCD 
phase boundary, location of critical point (CP) and determine the order of 
the phase transition. 



2 EOS for QGP and HG 

The EOS for QGP in our quasiparticle framework and the calculations 
regarding thermodynamical quantities like pressure, energy density, particle 
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Figure 1: Variations of trace anomaly (e — 3p) /T^ with respect to tem- 
perature at /xs = and 500 MeV in our hybrid model. Shaded portion 
represents the calculation of Andronic et al. J,. 



density etc. can be found in our earlier work [3l|4]. In this model, we start 
with the definition of average energy density and average number density 
of particles and derive all other thermodynamical quantities from them in 
a consistent manner. The parameters used in the calculations are given 
there. 

Recently we proposed a new excluded-volume model for the hot and 
dense HG [Sj O |6] . The grand canonical partition function in our excluded 
volume model of HG can be explicitly given as follows : 



Btt^T Jyo Jo yFT^ [exp (^^j + 1] 

where gi is the degeneracy factor of ith species of baryons,i?i is the en- 
ergy of the particle {Ei — y^k'^ + mf), V'^ is the eigenvolume assigned to 
each baryon of ith species and hence ^i^j becomes the total occupied 
volume where Nj represent the total number of jth baryons. 



3 Results and Comparison with Lattice QCD 

In Fig. 1, we have plotted the results obtained for the trace anomaly 
factor (e — 3p) jT^ in our hybrid model calculations using HG and QGP 
equations of state separately at /is = 0. We further compare our results 
with the results obtained from a recent lattice calculation We notice 
that our results yield an excellent fit to the lattice data. In Fig. 2, we 
show the variations of baryonic succeptibility normalized as x^/^^ with 
temperature at [Ib = and compare our results with the lattice QCD 
results. Our hybrid model results again compare well with the lattice data 
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Figure 2: Variation of normalized baryon number succeptibility with re- 
spect to temperature at fj,B = in our hybrid model. Lattice data points 
at /Lts = are taken from Ref. [5]. 
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Figure 3: Variation of speed of sound with the temperature at /is = MeV 
as obtained in our hybrid model and compared with lattice results [9l [10] . 
Dash-tripple dotted curve shows the result of Andronic et al. [7] . 




Figure 4: Variation of (As/T^) = [s/T^)qgp - {s/T^)hg with respect to 
coordinates of various phase transition points on the {T^^b) phase bound- 
ary. We have used transition points from Ref. [4]. 

points. Here again the curves, for 'X^ jT^ obtained for HG and QGP phases 
are smoothly connected at around T = 170 MeV, which shows the presence 
of a cross-over transition between two phases. In Fig. 3, we show the 
variations of square of speed of sound (c^) with respect to temperature at 
/iB — MeV and again a comparison is given with the lattice QCD results. 
For fiB ~ 0, our results reproduces the lattice results very well. We have 
also shown separately the curve for obtained by Andronic et. al. by 
using RGSG model for HG phase only [7 . However, the features of the 
curve differ from the lattice data, although it also yields a minimum at 
around the same temperature. 

4 Critical Point (CP) and Order of Phase 
Transition 

In this section we attempt to show the precise location and nature of CP 
existing on this phase boundary. 

In Fig. 4, we show what will happen to a quantity depicting the change 
in the entropy density from the phase transition at CP of the phase diagram. 
We define the normalized difference ^ = {s/T^)qgp — {s/T^)hg and 
demonstrate its variations with respect to the coordinates of the phase 
transition points lying at the deconfining phase boundary. We find that 
^ ^ 0.0 and positive along the deconfining phase boundary in the case 
of first order phase transition which supports the role of the presence of a 
nonvanishing latent heat in the phase transition from HG to QGP. However, 



we surprisingly notice that ^ « exactly at the CP and thus CP can be 
taken as a point where the first order phase boundary terminates and phase 
transition changes its order. 



5 Summary 

Thus above results give a firm indication that the order of phase transition 
changes at CP for the deconfining phase transition. We hope that our re- 
sults will clarify the mist surrounding the understanding of the deconfining 
phase transition and the conjectured phase boundary. More importantly, 
we have formulated a phenomenological hybrid model which provides a re- 
alistic EOS for the entire QCD matter and in the absence of a first-principle 
lattice QCD calculation especially at finite /is, it can be reliably used for 
deriving the information on the QCD phase boundary. [1] [2] . 
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